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Abstract. The explicit calculation of the scaling form of the two-time autocorrelation 
function in phase-ordering kinetics and in those cases of non-equilibrium critical 
dynamics where the dynamical exponent z — 2 through the extension of dynamical 
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extension from the usually considered d-dimensional ageing or Schrodinger algebras to 
a new kind of representation of the conformal algebra in d + 2 dimensions. Explicit 
tests in several exactly solved models of simple magnets and through simulations in 
the 2D Ising and g-states Potts models (q = 2, 3, 8) quenched to T < T c are presented 
and the extension to systems with non-equilibrium steady-states is discussed through 
two exactly solvable models as well. In the context of surface growth models, possible 
generalizations for a dynamical exponent z = 4 and beyond are discussed. 
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1. Introduction 



Phase-ordering kinetics [T7j may occur if a magnet is rapidly brought ('quenched') 
from some initial disordered state to some temperature T < T c below its critical 
temperature T c > 0. At equilibrium, such systems have at least two equivalent 
macroscopically ordered thermodynamic states. Because of the competition between 
these two equilibrium states, the system cannot simply relax within a finite time towards 
one of them but rather, stimulated by the small effective magnetic fields which are 
created by the environment of each magnetic atom, ordered domains will form and 
subsequently grow. Although the average magnetization, if initially zero, will not 
change, there is a slow, non-exponential relaxation as the walls separating the ordered 
domains slowly move. For a spatially infinite system, this gradual coarsening process will 
continue indefinitely, and the state of the system and the behaviour of time-dependent 
observables will depend on the age of the system, that is the time elapsed since the 
quench. In addition, as will be discussed in detail in this article, dynamical scaling 
is observed, although neither of the equilibrium states by itself is scale-invariant. The 
three properties (i) slow, non-exponential dynamics, (ii) breaking of time-translation 
invariance and fiii) some kind of dynamical scaling are the constitutive properties of 
ageing systems]^ 

Physical ageing was first identified in celebrated experiments on the mechanical 
response of certain glassy systems to external stress [97]. For a better theoretical 
understanding, the quite recent realization that very similar phenomena can also be 
found in simple magnets, without disorder nor frustrations, might lead to conceptual 
insights which in turn could become also fruitful in more complex systems. For recent 
reviews of this intensively studied topic, see HSl EDI I2H I2SI SD ESI ED E2] - 

To begin, we recall some facts about dynamical scaling of phase-ordering systems. 
Unless explicitly stated otherwise, we assume that throughout the order-parameter is 
non-conserved by the dynamics and that the initial state is totally disordered. For late 
time, linear size L of the ordered domains grows as L ~ t l / z and it is known that the 
dynamical exponent z = 2 for a non-conserved order-parameter [18]. If 4>(t,r) is the 
order-parameter at time t and location r the following dynamical scaling of the two-time 
autocorrelation C and (linear) autoresponse function R is usually assumed 

C(t, s) = <0(t, r)cf>(s, r)> ~ s- b f c (t/s) (1) 

~ ^- a f R {t/s) (2) 

h=0 

where h(s, r) is the conjugate magnetic field at time s and location r. The scaling 
behaviour is expected to apply in the so-called ageing regime where 



R(t, S )- '<«'■'» 



8h(s, r) 



t, S 3> tmicro and t S 3> t m icro (3) 

f In distinction to chemical/biological ageing, 'physical ageing' as defined here can be observed even 
if the underlying microscopic dynamics is completely reversible. 
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where t m i cro is a microscopic time-scale |lU4j . In writing eqs. (11121) it was tacitly assumed 
that the scaling derives from the algebraic time-dependence of the single characteristic 
length-scale L = Lit) ~ t 1 ^ which measures the linear size of correlated or ordered 
clusters and where z is the dynamic exponent. Under these assumptions, the above 
forms define the non-equilibrium exponents a and b and the scaling functions fc{y) and 
fniy). For large arguments y = t/s — > oo, one generically expects 

My) ~ y' Xc/z , f R (y) ~ y~ Xn/z (4) 



where Ac and Xr, respectively, are known as autocorrelation [391 EE] and autoresponse 
exponents [86]. In non-disordered magnets with short-ranged initial conditions one 
usually has Ac = Xr, but this is need no longer be true if either of these conditions 
is relaxed [HH [861 EU EE]. Field-theoretical considerations for systems quenched to 
T = T c show that for a non-conserved order-parameter the calculation of Xc,r requires 
an independent renormalization and hence one cannot expect to find a scaling relation 
between these and equilibrium exponents (including z) [68] |j] For phase-ordering 
(quenches to T < T c ) a well-known result due to Bray [T7] relates the difference Ac — Xr 
to an exponent describing the fall-off of the initial correlator. In particular it can be 
shown that Ac = Xr for short-ranged initial correlations [17] . This same result can 
alternatively be derived from local scale- invariance [87] (see section 3). 

The values of the exponents a and b are as follows for phase-ordering. First, one 
usually observes simple scaling of C(t, s) = fc(t/s), hence b = oJ§| Second, the value of 
a depends on whether the equilibrium correlator decays exponentially ('short-ranged') 
or algebraically ('long-ranged'). This defines the classes S and L, respectively and one 
has [El E6l EH El] 

J e-l r l/ f J class S J l/z ( 

° eqlrj ~ | | r |-( d - 2 +") ^ \classL ^ ° " \ (d - 2 + rj)/z {b) 

Examples for short-ranged models (class S) include the Ising or Potts models in d > 1 
dimensions, while the spherical model or the 2D XY model below the Kosterlitz- 
Thouless transition are examples for long-ranged systems (class L)jjJ 

A much harder problem is posed by the quest to find the universal scaling 
functions and fc(y)- Indeed, one might consider the analogy with equilibrium 

critical phenomena where simple scale-invariance can be extended to conformal 
invariance, under quite general conditions. Since conformal transformations are 
scale-transformations with a spatially varying rescaling factor (such that angles are 
conserved), one can inquire about the possibility to similarly extend dynamical scaling 
to an invariance of the physical system under more general, local scale-transformations 



| That is different if a non- vanishing initial value of the order-parameter is considered 
§ Explicit results on the ageing behaviour of the spherical model close to free surface give b ^ [7J. 
|| For quenches to T — T c , all systems are of course in class L. 



Autocorrelation functions 



4 



[4"9"1 I5T)] . Specifically, motivated by the analogy with conformal invariance [SHI [13 , the 
transformations of time 



were considered to which transformations of space r i— > r' must be added in such a way 
that, say, these local scale-transformation form a Lie algebra. Furthermore, it can be 
argued that the linear response functions should transform covariantly under such local 
scale-transformations (at least if the scaling operators from which they are built are so- 
called quasi-primary operators) [IHl EH [87] . Taking into account that time-translations 
cannot be part of local scale-transformations as applied to ageing phenomena (hence 
(3 = in (ED), this leads to [871 E7J [52] (Q is the Heaviside function) 



such that the form of fniy) is determined by the two exponents \r/z and a' which are 
related to the scaling dimensions of the order-parameter and the associate response 
field <p and fo is a normalization constant. By now there exists quite a long list of models 
where either the exact solution or numerical data for R(t, s) = (<p(t)<f>(s)) are compatible 
with the prediction (jTJ of local scale-invariance as reviewed in detail in [61, 62]. 

Here, we shall address the question how to find fc(y) from an extension of 
dynamical scaling to local scale-invariance. Although the method has been successfully 
applied in various contexts, see [551 El EH ETJ [57J [901 E], only the barest outline of it 
was ever published [551 E7]. As we shall see, the analytical calculation of correlators is 
quite a demanding task, be it in the context of closed- form approximations [T6l 11001 [91]. 
perturbative schemes [7H [7_E1 [79] or else from a dynamical symmetry. Since z = 2 in 
phase-ordering kinetics, a natural candidate for a local scale-symmetry is provided by 
the well-known Schrodinger group and/or some of its subgroups. We shall recall their 
definition in section 2 and remind that the Schrodinger group is a group of dynamical 
symmetries of the free diffusion equation. It may be less well-known that, provided one 
considers the diffusion constant as an additional variable, there is a larger dynamical 
symmetry isomorphic to a conformal group in d + 2 dimensions [201 E2] . These results 
may be extended to slightly more general diffusion equations. Next, we consider the 
Langevin description for phase-ordering kinetics where there exists an exact reduction 
formula which reduces the calculation of any n-point correlation/response function to 
the calculation of certain correlators/responses in the noiseless part of the theory only, 
given that so-called Bargman superselection rules for that deterministic part are valid. In 
section 3, we recall first the standard result from Schrodinger-invariance for the required 
three-point function and then show how using the previously discussed extension to a 
conformal group the two-time correlation function can be found. In section 4 we review 
existing tests of this prediction, both analytical and numerical. Finally, in section 5 we 
present a short outlook how the methods reviewed here might be extendable to cases 



* i— »> r = 



at + /3 
it + 5 



ad — /?7 = 1 



(6) 




(7) 
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where z 7^ 2 and conclude in section 6. A technical point on the determination of the 
autocorrelation function is treated in the appendix. 

2. Dynamical symmetries: Schrodinger, ageing and conformal 

In this section, we recall those elements of local scale-invariance, specialized to the case 
z = 2, which we shall need for the explicit calculation of two-time correlations. 

2.1. Schrddinger-invariance of the free diffusion equation 

Consider the dynamical symmetries of the free diffusion (or free Schrodinger) equation 

2Md t (j) = A(f) (8) 

where A = V • V is the spatial laplacian and the 'mass' M. plays the role of a kinetic 
coefficient. The Schrodinger- group Sch(d) (actually, already found by Lie more than a 
century ago) contains the space-time transformations 

. at + (3 , Tlr + vt + a 
t^t = ; r ^ r = — — , a8 - ^7 = 1 9 

jt + jt + 

where a, (3, 7, 5, v, a are real (vector) parameters and 1Z is a rotation matrix in d spatial 
dimensions. The group acts projectively on a solution of the diffusion equation through 
(t,r) h+ g(t,r), <}>^T g (f> 

(T g( f>) (t, r) = f g {g-% r)) <f>{g~\t, r)) (10) 

where g is an element of the Schrodinger group and the companion function reads [82|[85] 

M 7T* 2 + 2TZr ■ (7a - 5v) + 7a 2 - tSv 2 + 27a ■ v 



f g {t,r) = { 1 t + 5)- d l 2 



cxp 
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It is then natural to include also arbitrary phase-shifts of the wave function within 
the Schrodinger group Sch(d). We shall denote the Lie algebra of Sch(d) by sci) d . The 
invariance of the space of solutions of the free diffusion equation (|S]) under this group 
can be illustrated in d — 1 by introducing the Schrodinger operator 

S := 2M y_ 1 - Y 2 1/2 = 2Md t - d 2 r (12) 

In what follows, we shall restrict calculations to the case d — 1 since the extensions to 
d > 1 will be obvious. The Schrodinger Lie algebra sd) 1 = (X-x j o,x,Y_i i,M Q ) 
is spanned by the infinitesimal generators of temporal and spatial translations 
(X-i,Y-i/2), Galilei-transformations (Y72), phase shifts (Mo), space-time dilatations 
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with z = 2 (X ) and so-called special transformations (Xx). The generators read 
explicitly [19] 

x n = - t n+ % - ^t n rd r - ^±±W-v - |( n + i)* n 

F m = - r +1/2 9 r - (m + A<t m - 1/2 r (13) 

m„ = - 

Here x is the scaling dimension and At is the mass of the scaling operator on which 
these generators act. The non-vanishing commutation relations are 

[-^nj-^n'] = {n — n')X n+n i , [X n , Y m ] = ^— — rnj Y n+m 

[X n , M n ] = - n'M n+n i , [Y m , Y m >] = (m - m')M m+m < (14) 

The invariance of the ID diffusion equation under the action of scl) 1 is now seen from 
the following commutators which follow from the explicit form (flBl 

[S.ALJ = [<S,y ±1/2 ] = [5, M ] = 

[S,X ] = -S , [«S,Xi] = -2*5 - (2x - 1)M (15) 

Therefore, in d spatial dimensions, /or any solution <fi of the Schrddinger equation 
S(f) = with scaling dimension x = d/2, the infinitesimally transformed solution X(j) 
with X G sch d also satisfies the Schrddinger equation SXcf) = [70l |82| US]. 

Technically, we have been making an important assumption here. Indeed, and 
borrowing terminology from conformal invariance |T2], only so-called quasiprimary 
scaling operators 4> will transform under the action of the Schrddinger group such that 
the change under an infinitesimal transformation, generated by X e sd) d , is simply given 
by X(p jl9l ED] • It is non-trivial that a given physical observable should be represented 
by a quasiprimary scaling operator, although this appears usually to be the case for the 
order-parameter. For example, if <fi = <p(t, r) is quasiprimary, neither of the derivatives 
dt<f) nor d r <j) is. 

Now the form of n-point functions built from quasiprimary scaling operators can be 
constrained by writing down differential equations X^((f)\ . . . 4> n ) = with the n-particle 
extensions X™ of the generators X E 5ch d . Explicit results will be quoted below. 

2.2. Ageing invariance 

For applications to ageing, we must consider to so-called ageing algebra 

a Q t 1 := (X ,u Y -^ M o) C scfh (16) 

(without time-translations) which is a true subalgebra of sd) 1 . The generators Y m ,M n 
retain their form, but the generators X n only exist for n > 0. They now read [87J [57] 

X n = -t n+l d t - ^±t n rd r - (ra + 1)n At*"-V 2 - + 1)*" - int n ; n > (17) 

Zi 4r Zi 
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where £ is a new quantum number associated with the field <fi on which the generators 
X n act. When £ 7^ 0, the generators X n are only part of a representation of 
agc 1 but not of scf^. This is only possible for systems out of a stationary state 
(otherwise, the requirement of time-translation invariance and = 2X would 

lead to £ = 0). With respect to (TI5]) . only one commutator changes and we have 
[S, Xi] = —2tS — (2x + 2£ — 1)Mq, resulting in the weaker condition x + £ = 1/2 in 
order that the representation of agc 1 can act as a dynamical symmetry. It can be shown 
that in the scaling form eqs. ( 121171) of the two-time autoresponse function the difference 
a' — a is related to £ (and £ of the response field 0) and vanishes if £ + £ = 0. 

The meaning of this new quantum number £ of a quasiprimary scaling operator 
becomes apparent when we consider the extension [4~9l 157] of agt 1 to the infinite- 
dimensional algebra ao spanned by (X n ) ngNo , (Y m ) meZ+ y 2 and (M n ) neZ , such that 
eq. (tHJ) remains valid. It is well-known that Y m and M n generate time-dependent 
translations (which an additional phase) and phase shifts, respectively. On the 
other hand, the generators (X n ) n > eq. (jTTJ) are the infinitesimal generators of the 
transformation (t, r) i— > (£', r') where 



t = m , r = r<» (« 



such that /3(0) = and the scaling operator transforms as 



,'d(3{t')Y x/2 /dln/?(t')V* 
^• r )=l^ i UlnV 6XP 



A^r' 2 d / d/3(t') 
In 



4 df V dt 



0'(t',r') (19) 



This transformation is different from the one of a Schrodinger quasiprimary operator if 
£ 7^ but it suggests to define the scaling operator 

$(t, r ) := t~ e 0(t, r) (20) 

which is indeed quasiprimary under the ageing algebra at) and even the analogous 
extension so of the Schrodinger algebra sci) 1 , but with a modified scaling dimension 
x$ = + 2^ [57]. 

Therefore, in order to describe the local scale-invariance of ageing systems with 
z — 2, it is enough to study first Schrodinger-quasiprimary operators $ and only at the 
end, one goes over to ageing scaling operators (ft using eq. (120]) . 



2. 3. Conformal invariance of the free diffusion equation 

Important information will come from considering the 'mass ' M. not as some constant, 
but rather as an additional variable. It is useful to work with the Fourier transform of 
the field and of the generators with respect to M. and to define a new field ip as follows 
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</>(t,r)=<f> M (t,r) = -^= I dCe- iM ^(C,t,r) (21) 



In 

Provided lim^j-oo t, r) = 0, the diffusion equation becomes 

2i lWc + l£ = (22) 

which for brevity we shall also call a diffusion/ Schrodinger equation. In d = 1 the 
generators read (an eventual generalization to age x or at> is straightforward) 



Xn = - t n+1 d t - ^±1 t n rd r -(n+ l)*t n + i n{n + l \ n - l r 2 d f 



2 ' v '2 



Y m = - t m+1 ' 2 d r + i (m + t m ^ 2 rd c 
M n = it n d ( (23) 



This change of variables trades the complicated phases acquired by the field under a 
Schrodinger /ageing transformation for a time- dependent translation of the new internal 
coordinate (. For the Galilei-transformation Y\ii this was first observed by Giulini |42j. 

The passage to the new variable ( implies an extension of the dynamical symmetry 
algebra [52] which we illustrate for d = 1. Rewrite the physical coordinates (, t, r as the 
components of a three-dimensional vector £ = (£-i,£cb£i) where 

i=~(-£o + i£-i) , C = ^o + i£-i) , r = y|& (24) 

and i[)(£,t,r) = Then the diffusion equation (j22p becomes a three-dimensional 

massless Klein-Gordon equation 

3 M 3"tt(£) = (25) 

The Lie algebra of the maximal kinematic group of this equation is the conformal algebra 
conf 3 , with generators 

P =8 

= £pd v - Zvdp 

K/j, =2^ u d„-U v d lx + 2x^ (26) 
D = i v d v + x 

(/jl, v G {—1,0,1}) which represent, respectively, translations, rotations, special 
transformations and the dilatation. Hence the generators of sci) 1 are linear combinations 
(with complex coefficients) of the above conf 3 generators. Hence, for any d > 1, one has 
an inclusion of the complexified Lie algebra (scfj d )c into (conf d+2 )c [20J. Explicitly, for 
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A l 



W V + Xj 



V NX Y 

_ 1/2 



"®- 



X l . Y l/2 _ M 



T 



Figure 1. Roots of the complex simple Lie algebra B2 and their relation with the 
generators of the Schrodinger algebra sct) 1 and the conformal algebra conf 3 . The double 
circle in the centre denotes the Cartan subalgebra (). 



X_i = i(P_! - iP ) , X = -Id + ^M_ 10 , X, = -\{K^ + iK ) 

^-1/2 = -^? Pi , yi=-y^"(M_n + iMoi) , M = P_i + iP . (27) 

The four remaining generators needed to get the full conformal Lie algebra (conf 3 ) c can 
be taken in the form 

iV = —tdt + C^c time-phase symmetry 

VL = — (d r + ird t 'dual' Galilei transformation 

W = — ( 2 d^ — (rd r + i^dt — x( 'dual' special transformation 

V + = —2trd t — 2(rd^ — (r 2 + 2i(t)d r — 2xr transversal inversion. (28) 

The generators V_ and W are, up to constant coefficients, the complex conjugates of Yy 2 
and Xi, respectively, in the coordinates hence their names. The complex conjugation 
becomes the exchange t <-> ( in the physical coordinates ((, t, r). 

These results are illustrated in the root diagramme figured] [52J . It can be shown [72] 
that the generators considered above are roots of the complex Lie algebra B 2 — (conf 3 )c- 
Then to each root X e B 2 one may associate a two-dimensional vector f e Z 2 and linking 
the origin to one of the points in the root diagramme such that forming the commutator 
[Xi, X 2 ] = 012,3^3 corresponds to vector addition ji + p 2 = £3 and 012,3 = if the result 
of that addition falls outside the root diagramme. 



Autocorrelation functions 



10 



Since possible isomorphisms are described by the Weyl group which is generated by 
the simple symmetries in figure UJiui : (ei, e 2 ) ( — e 2> ~ e i) an d w 2 : (ex, e 2 ) 1— > (— ei, e 2 ) 
one may readily write down the list of non-isomorphic maximal subalgebras of B2 . These 
are [52] 

(i) the algebra agt l := age 1 © CiV, obtained from agt l by adding the generator N. 
This algebra is known as the minimal standard parabolic subalgebra^ 

(ii) the extension sci) 1 := sd) 1 (BCN of the Schrodinger algebra which is also a parabolic 
subalgebra. 

(iii) the parabolic subalgebra algebra alii := olti © CA" where alti : = 
(D, Y_i/2, M , Y1/2, V + , Xi). See [60j for a study of its representation theory 

Against often-read statements, which go back to [4|, sd) l cannot be obtained through 
group contraction from taking a non-relativistic limit of the conformal algebra conf 3 . 
Rather, this limit yields a projection conf 3 — > atti [52]. On the other hand, alti and its 
infinite-dimensional extension can be obtained from a group contraction, realized as a 
non-relativistic limit, of a pair of commuting Virasoro algebras [60J. 

These new generators are indeed dynamical symmetries of the free diffu- 
sion/Schrodinger equation. To see this, consider the ID Schrodinger operator 

S = M X_! - Y 2 1/2 = 2id c d t + d 2 r (29) 

It is straightforward to check that 
[S,V-] =0 = [S,N] 
[S,V+] = 2(1 - 2x)d t - 4rS 

[S,W] = i(l - 2x)d r - 2(S (30) 

which complement eq. ( flBI) . Hence, any solution ip of the ID Schrodinger equation 
Sip = with a scaling dimensions x = 1/2 is mapped onto another solution of 
under the action of the conformal algebra (conf 3 )c- Of course, the space-time realization 
of this conformal symmetry is unconventional in the sense that space-time rotations with 
M. fixed do not occur. 

A similar extension can be made for higher spin fields. The simplest example treats 
the relationship of the non-relativistic Levy-Leblond equations to the massless Dirac 
equation [59] . 

For extensions to d > 1 and the emerging mathematical theory of central extensions 
of the Schrodinger- Virasoro algebra 50 and its vertex operators we refer to the litterature 
[891 [TOT]. 

Given the conformal invariance of the free Schrodinger equation, the form of the 
two- and three-point functions is well-known [88J. Transforming these results back to a 

\ A parabolic subalgebra is defined as the sum of the Cartan subalgebra and of the positive root spaces 
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fixed value of Ai, one recovers the two-time and three-time response functions, together 
with an explicit proof of their physically required causality properties |52j . 



2-4- Reduction from Langevin to deterministic equations 

So far, we have considered the dynamical symmetries of simple diffusion equations but 
we should really have been interested in the possible symmetries of stochastic Langevin 
equations. For a non-conserved order-parameter (model A) one usually writes [63] 

where V is the Ginzburg-Landau potential and rj is a gaussian noise which describes 
the coupling to an external heat-bath and the initial distribution of 0. At first sight, 
there appear to be no non-trivial symmetries (besides the obvious translation and 
rotation invariance and possibly dynamical scaling), since the noise r\ does break Galilei- 
invariance. This can be seen from an easy calculation, but in order to understand this 
physically, consider a magnet which is a rest with respect to a homogeneous heat-bath 
at temperature T . If the magnet is moved with a constant velocity with respect to the 
heat-bath, the effective temperature will now appear to be direction-dependent, and the 
heat-bath is no longer homogeneous. 

Still, there are hidden non-trivial symmetries which come about as follows [87]: split 
the Langevin equation into a 'deterministic' part with non-trivial dynamic symmetries 
and a 'noise' part and then show using these symmetries that all averages can be 
reduced exactly to averages within the deterministic, noiseless theory. Technically, 
one first constructs in the standard fashion (Janssen-de Dominicis procedure) [32j [69] 
the associated stochastic field-theory with action J[4>, 0] where is the response field 
associated to the order-parameter 0. Second, decompose the action into two parts 

J[cj)A} = J [4>A}+J b {$\ (32) 

where 

Jo[0, 4>}= I dtdr \2Md t <f> - A0 + j^j (33) 



-Xll 



contains the terms coming from the 'deterministic' part of the Langevin equation (V is 
the self-interacting 'potential') whereas 

J b $\ = -T I dtdr 0(t, r) 2 - - [ dr dr' 0(0, r)a(r - r')0(O, r') (34) 

</R + xR d 2 J R 2d 

contains the 'noise'-terms coming from ( 13T1) [69]. It was assumed here that (0(0, r)) = 
and a(r) denotes the initial two-point correlator 



a(r) := C(0, 0; r + r', r ) = (0(0, r + r')0(O, r')> = a(-r) (35) 
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where the last relation follows from spatial translation/rotation-invariance which we 
shall admit throughout. 

It is instructive to consider briefly the case of a free field, where V = 0. Variation 
of ( 132]) with respect to and 0, respectively, then leads to the equations of motion 

2Md t (j) = A0 + T4> , -2Md t 4> = A0 (36) 

The first one of those might be viewed as a Langevin equation if is interpreted as 
a noise. Comparison of the two equations of motion (1361) shows that if the order- 
parameter is characterized by the 'mass' M. (which by physical convention is positive), 
then the associated response field is characterized by the negative mass —M.. This 
characterization remains valid beyond free fields. 

We now concentrate on actions J7"o[0, 0] which are Galilei-invariant. This means 
that if (.)o denotes the averages calculated only with the action j7o, the Bargman 
superselection rules [3] 

- 1 K,m (37) 

hold true. It follows that both response and correlation functions can be exactly 
expressed in terms of averages with respect to the deterministic part alone. The simplest 
example is the two-time response function (we suppress for notational simplicity the 
spatial coordinates) [HIj 

= (0(t)0( S )) = (0(t)0( S )e^) =(0(t)0(s)) (38) 



5h(s) 



where the 'noise' part of the action was included in the observable, the exponential 
expanded and the Bargman superselection rule fl37|) was used. In other words, the 
two-time response function does not depend explicitly on the 'noise' at all and this 
explains why the prediction (JTj) which was derived from the dynamical symmetries 
of a deterministic equation successfully describes the response functions calculated in 
stochastic models, as reviewed recently in [EH E2] . 

The correlation function is reduced similarly [87J 



\2 ' 

' o 



C{t, s;r)=T I dudR U(t, r + y)<f>(s, y)<f>(u, R) 

JR + xR d ^ 

+ 1 I dRdR' a{R- R') U{t,r + y)<f>{s,y)0{O, R)$(0,R') ) (39) 
2 J W 2d \ I o 

Only terms which depend explicitly on the 'noise' remain - recall the vanishing of the 
'noiseless' two-point function (0(t)0(s))o = because of the Bargman superselection 
rule. Eq. ( 1391) is the basis for our later explicit calculation of C(t, s). 
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2.5. Extensions 

For the later applications to phase-ordering, we need the following extension of fl3Tj) 

2M^ = ^-*M_ B(t)# + , (40) 

where the time- dependent potential v(t) arises naturally in spherical models. In more 
general models, it might be viewed as a Lagrange multiplier to enforce that C (t, t) is 
finite. Through the gauge transformation 

c(>(t, r) ^ c/>(t, r) k(t) ; k(t) := exp JL J\ t ' v(t')^ (41) 



eq. flU]) is reduced to (Ell) [HZ]. 

Finally, we consider semi-linear equations ( |3T|) where V = V[0, </>]. Here the 
essential ingredient is to recognise that the coupling constant g of a non-linear term 
will in general have some scaling dimension y g . This implies in turn that under local 
scale-transformations g must be transformed as well and consequently, the infinitesimal 
generators have to be generalized accordingly. Standard methods [15] then allow to 
compute the Schrodinger-invariant semi-linear equations and one has: 

(i) for fixed masses, the invariant semi-linear Schrodinger equations are of the form 

S<j> = <j ) U^) 1/X f(g x (^Y 3 ) (42) 



where x is the scaling dimension of <fi and / an arbitrary differentiable function [6]. 

(ii) for variable masses, the ID invariant semi-linear Schrodinger equations with real- 
valued solutions are of the form 

S^ = ^f{g^) (43) 

where / is an arbitrary differentiable function [96J. We point out that in this way 
one can obtain Schrodinger-invariant diffusion equations with real- valued solutions, 
at the expense that after inverse Laplace transformation with respect to ( products 
if) ■ if) ■ . . . are replaced by convolutions <$*<$*... with respect to A4. 

Slightly more general forms are possible if only ageing-invariance is required [SI ES] ■ 



3. Calculation of two-time autocorrelation functions 

3.1. Schrodinger-invariant autocorrelators 

We now turn to the explicit calculation of the two-time autocorrelation function, starting 
from the exact reduction formula eq. (1391 . As initial condition, we use a fully disordered 
initial state, with a(r) = aoS(r), see eq. ( |35l) . Then the calculation reduces to find 
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a three-point response function (0(t)0(s)0 2 (w)) o within the deterministic part of the 
theory and local scale-invariance will become useful here. At first sight, since there no 
time-translation-invariance, it might appear that the only dynamical symmetry available 
would be <XQt d , with rather limited predictive capacities. It is better to proceed as follows 
[551 E7j: 

(i) we break time-translation-invariance explicitly by considering a Langevin equation 
of the form (j4"0l) with a time-dependent potential v(t). We can reduce this to the 
standard form (I3TI) via the gauge transformation fj4~T]) 



k{t) ~ t F or equivalently v{t) = — ~J~ (44) 

(ii) recall that for ageing-invariance the scaling operators <fi (with scaling dimension x) 
do not transform as conventional quasi-primary scaling operators, but are related 
via fl20l) to bona fide quasiprimary scaling operators $, with scaling dimension x+2£. 

(iii) we then use full Schrodinger-invariance for the calculation of the required three- 
point function ($(t)$(s)S 2 (w)) . 

(iv) since this procedure still leaves a scaling function of a single variable undetermined, 
we use the further extension from Schrodinger to (d + 2)-dimensional conformal 
invariance. This requires to consider the masses as further variables, however. 

For illustration, we briefly consider the autoresponse function R(t, s) = (<f>(t)<f>(s))o. 
According to the procedure outlined above, <j) and <fi are characterized by the exponents 
x, £ and x, £, respectively, in addition to the parameter F coming from the gauge 
transformation (141144)) . This leads to [57J 



R(t,s) = (t^(t,r)s^(s,r 
k{t 



^ 2)/ MM (;-l) ®(t-s)5 x+2 ^ +2l (45) 



which reproduces (JTj), including its causality, with the following identification of 
exponents (we use z = 2 here) 

\ R = 2(x + £-F) , d-a = £ + £ , l + a = Ux + x) (46) 

Alternatively, the gauge exponent F may be expressed as F = \(x — x) + £ — £. 

Similarly, the autocorrelation function can be written as, where the composite field 
4> 2 is characterized by 2x 2 and 2£ 2 
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C(t,s) = ^ / dRR ( 3 \t,s,r ini] R) + ^- [ du [ dRRi 3 \t,s,u;R) (47) 



Ro ] (t, s, u; r) := ( <j>(t\ y)(/>(s; y)(j>{u\ r + y) 



(tsfu 2 ^ ($(t; y)$(s; r + y) 



i„2fa #( ^(3) 

k 2 {u 



(48) 



Here r; n i is some small 'initial' time which sets the beginning of the scaling regime. At 
a later stage, the limit r ini — > should be taken. 

Throughout, we consider $ 2 as a composite field with scaling dimension 2x 2 + 4^2 
(only for free fields, one would have x = x<i and £ = £2)- Schrodinger-invariance gives 
the three-point response function llf for v(t) = [49] 



7?.g 3 ' ) (t, s, u; r) = 7£q' 3 '' (£, s, w) exp 



(3), 



M t + s - 2u 
~2 (s~u)(t-u) 



2(t-u)(s 



u 



H { 3) (t, 8, u) = Q(t - u)Q(s -u)(t- u )-**-** (s - up 2 ' 26 (t - 3)-*-^+% (49) 

where \I> = ^/(p) is an arbitrary scaling function, to be determined. Of course, and 
we shall come back to this below, this is only applicable within the ageing regime as 
defined in eq. (131). It is useful to write down the autocorrelation function in the form 
C(t, s) = C t h(t, s) + C pr (t, s). Here the thermal part and the preparation part C pr 
(where the limit r; n i — ► was taken) are given by 



Cth(M) 



Co,th S 

X 



d/2+l-x-x 2 [ _ 
S 



1-1 

S 



X2+2&-x-2£-d/2 



t 

V 

s 



C pr (t, s) 



&(w) 



fi d/2+2F-x 2 -2&-x 



MW R 2 



d/2-x 2 -2^ 2 



dR exp 



l-i 

s 



X2-x+2&-2£-d/2 



t/s + 1- 2v 
t/s- 1 

t/s + 



t/s 



tf (K 2 ) 



(50) 



and we have explicitly used s < t. Here Co.th is a normalization constant proportional 
to the temperature T and Co iPr is a rescaled normalization depending on the initial 
time-scale r ini . 

For the interpretation of these forms, we concentrate on simple magnets relaxing 
towards equilibrium stationary states. From renormalization arguments [T?] |2T] one 
expects that the thermal part C t h will be leading one for quenches onto the critical 
temperature T = T c while the preparation part C pr should be dominant one for quenches 
into the ordered phase (T < T c ). Comparison with the expected scaling forms (fTf4"|) then 
give 
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(i) T = T c : for relaxation to equilibrium, one a = b and from (1501) we read off 

a = b = x + x 2 -l-^ , \ c = 2(x + £-F) (51) 
The autocorrelation function scaling function becomes, up to normalization 

C\y) = y 1+a '- Xc/2 (y-i) b - 2a ' +2 ^ 1 



(52) 



where /i := £ + £2 is a further parameter whose value is not predicted by the theory, 
(ii) T < T c : one has 6 = and from the asymptotic scaling function we find 

\ c = 2(x + Z-F) (53) 
The scaling function reads simply, up to normalization 

Since x 2 — x + 2£ 2 — 2£ = |<i — Ac < because of the Yeung-Rao-Desai inequality 
|103j . this exponent combination will only vanish for free field. This illustrates 
again that x 2 and £2 must be treated as independent quantities and cannot a priori 
be related to other exponents. 

Comparing the above results for Ac with eq. (H6|) giving A#, we see that for short- 
ranged initial conditions and a deterministic part J§ of the action which is Schrddinger- 
invariant, the autocorrelation and autoresponse exponents are equal, Ac = Xr [EI]- This 
well-known result [IT] can therefore also be understood as a consequence of a dynamical 
symmetry. 

Schrodinger-invariance by itself does not determine the form of ^(p). However, 
the related three-point response function should be non-singular as t — s — > 0. This 
requirement fixes the asymptotic behaviour for w — > 00 [871 [57] 

,tw \ -p ) X c ; if T < T c and 6 = 

' l + 2a>-b-2» ; if T = T c and a = 6 (55) 

This suggests the following approximate forms of the scaling function, with y = t/s 
f i(y+iy/(4y)]- Xc/2 ; if T < T c 

fc{y) « I yl+a'-Ac/Z ^^-2-2^ [(y _„)(!_ w )] a '^ 2 ^ ; if T = T c (56) 

[ x(y+ i-2v) b - 2a '- 1+2 ^ 

At least, these forms are consistent with the required asymptotic behaviour of fc{y) 
as y — > oo. For free-field theories (JS"5l) holds for all values of w and then ( ]56l) becomes 
exact. 



§ If one attempts to apply this kind or argument to systems relaxing to critical non-equilibrium steady 
states, one must take into account that a ^ b in general, see |61) and references therein. 
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3.2. Conformal invariance 

In order to find the last remaining scaling function St, we now propose to extend 
Schrodinger-invariance further towards conformal invariance in d + 2 dimensions as 
described in section 2. This implies that we must consider the 'mass' Ai of the scaling 
operator $ as a further variable. Consider three Schrodinger-quasiprimary scaling 
operators $ a = & a (t a , r a , -Ma) where a — a, b, c. Using the variables 



r = t a — t c , a = t b -t c , r = r a - r c , s = r b - r c 
the three-point function is 



($ a $ 6 $ c ) = 5 




F(t, a,r,s;M a ,M b ,M c ) 



(57) 



(5f 



where the delta function expresses the Bargman superselection rule and because of 
causality, r > 0, o > 0. Schrodinger-invariance almost fixes the form of F: 

M a r 2 M b s 21 



cr 



-73 



cxp 



a 



*(p,M a ,Mt 



(ra 



ST) 



and 



7i 



2ra(r — a) 



~%ac,b £,a 



72 



-Xbc ,a 



73 



(59) 
(60) 

(61) 



with x abjC := x a + x b - x c and ^ abyC :=£ a + £b-£c- 

The last function \I/ = ^ abc will now be found by extending invariance under sch d 
to invariance under the conformal algebra conf d+2 . For the explicit calculation, it is 
enough to consider the case d = 1. From the root diagramme figured] it is clear that 
if one has invariance under the generators iV and V_, then invariance under the other 
generators follows. The explicit single-particle form of the generators N and V- is, for 
d= 1 

N = -tdt-1- Md M 

V_ = - \d M d r + \rd t (62) 
Then the covariance of the three-point function under N and V- leads to 
(rd T + ad a + 2 + M a d Ma + M b d Mb ) F = 

(rd T + sd a - d Ma d r - d Ma d r ) F =0 (63) 

In turn, these imply the following conditions for the as yet unknown function \1/ 

[(2 - 7i - 72 - 7s) - pd P + M a d Ma + M b d Mb ]V = (64) 

" T S 

- [1 - 7i + MadM a + d Ma d P ] + - [1 - 72 + M b d Mb + d Mb d p ] 
It a 



+ 



r 



a 



[-73 - pd P - d Ma d p 



(65) 
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In eq. (I65p . it can be seen that if two of the inner square brackets vanish when applied 
to the vanishing of the third one follows. Consequently, we have the following linear 
system of three equations for the function \1/ which depends on three variables 

[(2 - 7i - 72 - 7a) - pd P + M a d Ma + M b d Mb ] V = 
[l-lx+M a d Ma + d Ma d p \y =0 (66) 

[1 - 72 + M b d Mb + d Mb d p ] =0 
The first of these leads to the scaling form 

KA KA\ 2-71-7-2-73 u ( M a + M b M a - M b \ 

ty{p,M a ,M b )=p 717 73 in p, pi (67) 

and the other two give the following system of equations for K = K(j], () 

[2 - 7l - 72 + (3 - 71 - 72 - 73)9, + r)d v + (d c + V d 2 + (d v d c ] K( V , () = 

[72 - 7i + (3 - 71 - 72 - 7 3 )<9 C + Vd< + (d v + + V d v d c ] K{ V , () = (68) 

which in principle can be solved through a double series in 7] and (. 

For our purpose, namely to find the two-time autocorrelation function, we merely 
need the special three-point function ($$$ 2 ) . In this case, we have 

7i = 72 = £2 + 26 , 73 = z-£ 2 + 2£-26 (69) 

M a = M b = M and M c = -2M. Hence * = = p 2 - x ^- 2 ^K{Mp, 0) 

where, up to normalization, K(rj, 0) is the solution of 

[2 - 2 7l + (3 - 2 7l - 73)9, + rjd, + rjd 2 ] K( V , 0) = (70) 

This is readily solved in terms of confluent hypergeometric functions and the final result 
for the scaling function \I/ reads, where ^0,1 are arbitrary constants 

*{p, M, M) = fa p 2 ~^-2?-26 iFi (2 _ 2 j 2 _ 4 ^ 3 _ x _ ~ _ 2 £ - 26; -Mp) 
+ fa M *+x2+^2~2 iFi ^ _ ~ + 2 ^ _ 2^ x + j 2 + 2 £ + 26 - 1; -Mp) (71) 

The main applications of eq. ( ITTi) will be to phase-ordering kinetics (T < T c ), where 
z = 2 for simple magnets without disorder [18] and furthermore 6 = 0. Then 

(p, M, M) = fa p x c~ 2a '- d /\F 1 (2A C - d - 2a', X c + l-2a'- d/2; -Mp) 
+ fa M 2a ' +d,2 - Xc l F 1 (X c - d/2, l + 2a' + d/2 - X c ; -Mp) (72) 

The case a = a! was already stated in [55]. On the other hand, for non-equilibrium 
critical dynamics (quenches to T = T c ) with dynamical exponent z = 2 and also a = b 
we have 

V(p,M,M) = fa p l - b ~ 2> *- dl2 1 F l (2 - d + 2a - 2b - 4p, 2 - b - 2p - d/2; -Mp) 
+ fa M^+^'hFt (1 - d/2 + 2a' - b - 2p, d/2 + b- 2/i; -Mp) (73) 
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We point out that the previouly derived result eq. (JTJ for the response functions does 
remain valid for this type of conformal invariance [52] and that the exponent a' of the 
autoresponse function plays a role in fixing the form of C(t, s). 

Before we insert these explicit results into the scaling forms eqs. ( 154"1) and ( 1521) . we 
should consider under which physical conditions their derivation is valid. In particular, 
we used from the outset dynamical scaling and and therefore restricted to the ageing 
regime eq. ([3]). Especially the condition t — s ^> t m i cro must be satisfied, see |104] for 
a careful discussion of this point. In our context, this means that for small arguments 
p — > 0, the form of the function \I/(p) is not given by local scale-invar iance. Rather, for 
p < 1 one should expect that the two-time autocorrelation function C(t, s) = C(s,t) 
should be symmetric in t and s and especially C(t, s) should be non-singular in the limit 
t — s — > 0. The requirement of the absence of a singularity in C(t, s) as t — s — > is 
equivalent to an analogous requirement on the three-point function ($(t)$(s)& 2 (u))o 
and leads to 

~ p?-3*+K-2b ; asp^O (74) 

Specifically, for the two physical applications under study, one recovers eq. ( 1551) . Since 
this condition is not satisfied by the solution given in eq. (1711) . one might consider the 
possibility that local scale-invariance only holds for sufficiently large arguments, say 
p > e and to write 

7 n ' y- (75) 

where ^lsi(p) is the expression given by eq. ( 17T1) and the constant e sets the scale which 
separates the two regimes. Finally, the constant \l/o is found from the condition that 
^(p) is continuous at p = e. 

In this way, one may obtain an explicit scaling function which at least satisfies 
the basic boundary conditions required for a physically sensible scaling function. For 
models with an underlying free-field theory, the approximation eq. (155]) becomes exact 
and leads to eq. (1561) . in agreement with the available exact solutions, as we shall review 
in section 4. For the case of phase-ordering, a straightforward, but just a little lengthy 
calculation (see the appendix for details) leads to 

= AT(\ C - 2a') 2 Fi(2\ c -d- 2a', \ c - 2a; X c +l-2a'- d/2; -I / w)w 2a ' ~ Xc 
+ BY (d/2) 2 Fi (A c - d/2, d/2; 1 + 2d + d/2 - \ c ; -l/w)w~ d/2 



+ AE 



l-2o' LA C 



2\ r -d- 2a' 



(Ew^'-^iXc + 1 - 2a', Ew) - 7 (A C , Ew) 



w- Xc 



X c + 1 - 2a' - d/2 

+ AE- 2a 'w- Xc ~i(X c , Ew) - A^iXc - 2d, Ew)w 2a '- Xc 

+ BE d/2 - Xc w- Xc -f(X c , Ew) - Bj(d/2, Ew) W - d/2 (76) 

where 2-F1 is a hypergeometric function and 7(0,2) an incomplete gamma-function [lj, 
E = M.E and A, B are constants related to ^ ,i- The first two lines would give the 
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z 


fc(y) 


ID Glauber-Ising, T = 


2 


(2/7r)arctan(v/2/(y-l)) 


spherical model, T <T C 


2 


((y + l) 2 /4y)- d / 4 


spherical, T = T c 


2<d<4 


2 


(y-iH+y-^ + i)- 1 


spherical, T = T c 


d>A 


2 


1 1 

( y _ l)-f +1 _ ( y + l)-f+l 


contact process bcp = EWl 


2 




pair 
contact 
process 


a < ac 


d>2 


2 


( y _i)-i+i_( y + i)-i+i 


a = ac 


2<d<A 


2 


( y + l)-| 2J p 1 (f ) f;f + l;^ T ) 


d > 4 


2 


(y + l)-§+ 2 - (y - l)-i+ 2 + (d- 4)(y - 


Edwards- Wilkinson ew2 


2 


( y _ i)-|+i+p _ ( y + i)-|+i+p 


Mullins-Herring MHl 


4 


(y _ (y + i)-f+i 


Mullins-Herring MH2 


4 


( y _ l)-|+l+p/2 _ ( y + l)-|+l+p/2 


Mullins-Herring MHc 


4 


( y _ l)-(*-2)/4 _ ( y + 1 )-(d-2)/4 



Table 1. Dynamical exponents and scaling functions (up to normalization) of the 
autocorrelator of some analytically solved models, see text for the definitions. 



scaling function if the above-described procedure were unnecessary. For free fields, 
one has Ac = d/2. Setting A = 0, we then have ty(w) = BT(d/2)w~ d ^ 2 and recover the 
case T < T c of fl55|) . 

Together, eqs. (1541176 j) give the final autocorrelation function C(t, s) of phase- 
ordering kinetics. This is the central result of this section (the special case a = a' 
was announced in |55j). 

We shall not give here the analogous expression for quenches to T = T c , since the 
known exactly solvable cases with z = 2 (such as the spherical model or the XY model 
in spin-wave approximation) are described by a free-field theory and in more general 
cases the dynamical exponent z ^ 2 (21] which in full generality is beyond the scope of 
this article. 

4. Tests 

We now compare the predictions of local scale-invariance - either eq. (156|) for models 
with an underlying free-field theory or else eqs. (I54II76I) or ( 15211761) if that is not the case 
- to explicit model results. In tabled] we collect results for fc(y) of some analytically 
solved models which we now proceed to analyse. 



Autocorrelation functions 



21 



4-1. One- dimensional Glauber-Ising model 

Our first example is the ID Ising model with Glauber dynamics [13] , at zero temperature 
and with nearest-neighbour hamiltonian H = —JY2i a i (T i+ii °i = =tl- F° r a totally 
disordered initial state, the exact two-time autocorrelation and autoresponse functions 
are 



C(M) = ^arctan^^- T , fl( t , a ) = ly (77) 

These results remain valid for long-ranged initial conditions (oyoo) ~ r~ v with 1/ > 
PI- 

Comparing _R(t, s) with the LSI-prediction eqs. (12171) one reads off a = 0, a' = —1/2 
and \r/ z = 1/2 [87]. Since the autocorrelation function does not depend on the range 
of the initial condition, it is reasonable to consider T = as the critical point of the ID 
Ising model. Since b = a = and Ac = \r = 1 and the ID Glauber model is described 
in terms of free fermions, we may use ( 155|) at T = T c = and then have from eq. ( |56j) 

pi V/f \ -1-2/1-1/2 /, \ 2/t 

C{t,s) = C j dvv 2 ^ (--l){l-v) U- + l-2v) (78) 

Choosing /x = —1/4 and C = V^/7r, the exact result quoted above is reproduced [57] . 



XY model in spin-wave approximation 

Consider a system with planar spins S G M 2 and with |S| = 1. In the spin- wave 
approximation, the usual nearest-neighbour hamiltonian becomes 

H = - S ( r ) ■ S ( r ') - \ [ dr WW)' (79) 

(rP) 2 jRd 

Since the spin-wave approximation considers the small deviations from a fully ordered 
state, LSI as formulated here with a vanishing initial order-parameter is not applicable. 
Therefore, a comparison of the existing results [M] for the spin correlators and responses 
with LSI would require the inclusion of a non-vanishing initial magnetization into the 
latter. On the other hand, since angular variable satisfies a simple noisy diffusion 
equation and has a vanishing initial value. The exactly known angular autoresponse 
and autocorrelation functions do agree, respectively, with eqs. ( I217I) and eq. (1561) in the 
critical case T = T c , where a — b — d/2 — 1 and Ac = Xr = d (for d > 2) [87\\]j 

|| If d = 2, A_r = 2 and A c = 0. 
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(7 
H 


T/T, 


Xr 


A 


B 


E 


undate 


Ref. 


2 





1.25(1) 


-0.601 


3.94 


0.517 


R 


m\ 




0.6610 


1.25(1) 


-5.41 


18.4 


1.24 


R 


m 




0.6610 


1.24(2) 


-5.41 


18.4 


1.24 


C 


m 


3 


0.5000 


1.19(3) 


-0.05 


2.15 


0.6 


C 


[76] 


8 


0.5001 


1.25(1) 


-0.07 


1.98 


0.4 


C 


[76] 



Table 2. Parameters describing the autocorrelation function of 2D Ising and g-states 
Potts models, according to LSI, eq. (f76|) . The letters R and C refer to random and 
checkerboard update, respectively. 

4-3. Spherical model 

The spherical model is described by a set of continuous spin variables (f) G R subject to 
the spherical constraint ^ r {4>{t, r) 2 ) = Af where M is the number of sites of the lattice. 
In the continuum limit, the equation of motion is [15] 

d t (j)(t, r) = Acj)(t, r) - v(t)<p(t, r) + V (t, r) (80) 

where 77 is the usual gaussian noise with zero mean and variance r)(t,r)rj(t,r')) = 
2T5(t — t')5(r — r'). The time-dependent potential v(t) arises from the spherical 
constraint. It can be found from the function git) = exp(2 J^drv (r)) where g(t) solves 
the Volterra integral equation, for short-ranged initial correlations 

g(t) = f(t) + 2T fdt' f(t - t')g{t') , /(*) = 0(t) (e- 4t I (U)) d (81) 

and Iq is a modified Bessel function [1] . This is exactly the kind of equations of motion 
of the type (140!) considered in section 3. Indeed, an explicit analysis shows [13 [2j [81] 
that g(t) ~ t r , with F = d/4 if T < T c . If T = T c , then F = 1 - d/4 for 2 < d < 4 
and F = if d > 4. It is therefore clear that the predictions ([7]) and fl56l) of local 
scale-invariance will be satisfied both below and at criticality (since z = 2 throughout) 
as indeed has been checked explicitly [87], see also table [TJ 

4-4- Two-dimensional Ising and Potts models 

All the exactly solved models reviewed so far are special in that the response and 
correlation functions satisfy linear equations of motion and therefore only test the rather 
weak prediction (1561) . It is therefore significant that there is evidence in favour of 
LSI (including the non-conventional conformal invariance described in section 3.2) from 
models whose equations of motion are non-linear. We refer to [EH [62] for the response 
function and concentrate on the correlation function for T < T c . 

To this end, the two-time autocorrelation function was calculated in the 2D kinetic 




Figure 2. Two-time autocorrelation function C(t,s) for the 2D Ising model and 
comparison with the LSI prediction (|54I76|) . In (a) the scaling behaviour for T = 
is shown for a large range of values of y — t/s and (b) redisplays the same data for 
smaller values of y. The inset compares the data for y w 1 and s = [200, 400, 800, 1600] 
with the closed approximation of BPT 16, 100J. The dash-dotted line in (a) is the 
second-order calculation from [75] and the dash-dotted line in (b) gives the result from 
[74], see text. The autocorrelation is re-plotted as a function of t — s and compared to 
LSI for T = in (c) and for T = 1.5 in (d). [55]. 



Ising model with heat-bath dynamics [55J as well as for the g-states Potts model, with 
q = 2, 3, 8 [7H1 EH [75]. The two-time autocorrelation was measured as 

C(t, S ) Ising = 1 , C(t, s) Potts = -L. ^ £<^ (t)i(TiW ) - 1 j (82) 

where the sums run over the entire lattice with M sites. In the first case, the <7j = ±1 are 
usual Ising spins, whereas in the second case (jj = 1,2, ... ,q are Potts spins. Another 
difference between the two simulations is that either random or checkerboard update 
was used. The results for the parameters describing C(t, s) from the two simulations are 
shown in tableland allow some explicit tests of universality. Throughout, the relation 
a = a' was assumed. First, the estimates for the exponent Ac in the q = 2 universality 
class, in agreement with earlier results [39], [19] , appear to be independent of the lattice 
realization and of the temperature, as long T < T c , which is of course expected from 
renormalization arguments [T7J. It is also satisfying to see that the same values of the 
parameters A, B, E of the LSI formula fl54|76j) can be used for two different update 
schemes. One observes that the term parametrized by B is apparently much more 
important than the other one, but the ratio A/B apparently depends on temperature. 

In figures [2] and [3] the degree of agreement between the numerical data and the 
prediction (1541176 j) is shown. Clearly, the system ages, as can be seen in figure [2bd, and 
there is excellent dynamical scaling for all models considered. For the Ising model, we 
show in figure E^ib the free-field approximation ([56]) . labelled 'app', which only describes 
the data for rather large values of y = t/s. This is hardly surprising, since the model is 
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Figure 3. Two-time autocorrelation function for the 2D g-states Potts model with 
q = 2 (left), q = 3 (middle) and q = 8 (right) [76l [67]- The different data symbols 
correspond to different waiting times. The solid lines show the fits to the scaling 
prediction (|76[) based on local scale invariance. 



known not to be described by a free field. On the other hand, we see that in all models 
the LSI predictions eqs. (I54II76I) give a good overall description over the whole range of 
the scaling variable, with the only exception of the relatively small region t/s < 2 — 3. 
It is conceivable that the patching procedure outlined in section 3 is not precise enough. 
It remains an open problem to derive a quantitatively precise formula which reproduces 
the numerical data over the entire range 1 < y < oo. 

On the other hand, the difficulty of achieving any quantitative agreement of an 
analytical approach and the numerical data is illustrated by comparing with the results 
of earlier attempts to find C(t,s) in phase-ordering kinetics. Using Gaussian closure 
procedures, Bray, Puri and Toyoki (BPT) find for the 0(n)-model [TBI 11001 [91] 

n [ R A n + lY| a / 4y \ d/i (l 1 n + 2 ( Ay \ d/2 \ 

where B is Euler's beta function [TJ . While this form does follow the general trend of the 
data, especially if y is not too large, but without reproducing them entirely, it mainly 
suffers from the disadvantage that it predicts the autocorrelation exponent Ac = d/2 
in disagreement with the data, see also figure [2b. Liu and Mazenko [73] and Mazenko 
[78] tried to remedy this by constructing more elaborate analytical schemes of which 
the BPT formula represents the lowest order. In particular, this leads for the 2D Ising 
model to fairly accurate values Ac = 1.246(20) [73] and Ac = 1.105. . . [78] (Mazenko 
also showed [79] that \r = Ac in his scheme). However, the detailed behaviour for fc{y) 
deviates largely from the data, as can be seen from the dash-dotted lines in figure [2b for 
the prediction of [73] - in agreement with the earlier tests performed by Brown et al. 
[T9] - and in figure [2h for the prediction of [7S] . These examples illustrate the extreme 
difficulty of describing C(t, s) theoretically. We consider it remarkable that local scale- 
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invariance achieves, for the first time, a precise representation of the numerical data 
over almost the entire range of the scaling variable y. 

A similar behaviour as for the Ising model is found for the Potts models with q = 3 
as well as q = 8 [76]. This is the first time that the ageing behaviour and LSI of the 
autocorrelator was tested for a system which undergoes a first-order transition. 

The examples presented here are the only ones which presently allow to test the 
extension a$z d — > sci) d — ► conf d+2 of the local dynamical symmetry. 



4-5. Bosonic contact and pair- contact processes 



Two exactly solvable models allow to test the of LSI in situations when the stationary 
states are no longer equilibrium states, which requires a modification of the reduction 
formulae to the deterministic part. See [61] and references therein for the motivation of 
these models, which are defined as follows [MJ E31 [5]. Consider a set of particles of a 
single species A which move on the sites of a <i-dimensional hyper-cubic lattice. On any 
site one may have an arbitrary (non-negative) number of particles. Single particles may 
hop to a nearest-neighbour site with unit rate and in addition, the following single-site 
creation and annihilation processes are admitted 

mA — (m + l)A , pA — —> (p — t)A ; with rates /i and A (83) 

where £ is a positive integer such that \£\ < p. Consider the following special cases: 

(i) critical bosonic contact process: (bcp) p = m = 1. Here only i = 1 is possible. 
Furthermore the creation and annihilation rates are set equal fi = A. 

(ii) critical bosonic pair-contact process: (bpcp) p = m = 2. We fix £ = 2, set 2A = ji 
and define the control parameter 

° - % < 84 > 

Linear equations of motion for one- and two-point functions may be derived 
straightforwardly [64j [HH [5] and it follows that in both models the spatial average 
of the local particle-density J dr p(t, r) = po remains constant in time. Scaling occurs 
along the critical line [83] 

£\ = /i. (85) 

The physical behaviour of these models is indicated in the phase-diagramme figure HI 
Whereas in the bcp the variance (p 2 {t, r)) diverges as t — ► oo for dimensions d < 2 but 

% For a coagulation process 2 A — > A, where I — 1, analogous results hold true if one sets A = \i and 
a = fi/D. 
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stays finite for d > 2 such that one has the same behaviour along the critical line, in 
the bpcp there is a critical value olq of the control parameter, given by 

i r°° 

— = 2 / du (e- Au UAu)) d (86) 

OLC Jo 

Specific values are «c(3) ~ 3.99, ac{&) ~ 6.45 and linid\ 2 oic{d) = 0. Now the variance 
(p 2 {t, r)) stays finite on the critical line for a < ac but diverges otherwise. 



X 





(b) 



Figure 4. Schematic phase-diagrammes for D ^ of (a) the bosonic contact process 
and the bosonic pair-contact process in d < 2 dimensions and (b) the bosonic pair- 
contact process in d > 2 dimensions. The absorbing region 1, where limt^oo p(x, t) = 0, 
is separated by the critical line eq. (|85|) from the active region 2, where p{x, t) — > oo as 
t — > oo. Clustering along the critical line is indicated in (a) and (b) by full lines, but in 
the bosonic pair-contact process with d > 2 the steady-state may also be homogeneous 
(broken line in (b)). These two regimes are separated by a multicritical point. 



This means that at the multicritical point at a = ac there occurs a clustering 
transition such that for a < « c the systems evolves towards a more or less homogeneous 
state while for a > ac particles accumulate on very few lattice sites while the other 
ones remain empty. In contrast with the bosonic contact process, clustering occurs in 
some region of the parameter space for all values of d. 

Solving the equations of motion with uncorrelated initial conditions yields the exact 
two-time autocorrelatoiE] (see tableland note that there is no scaling for a > ac) which 
may also be written in an integral form 



fc(y)=C I d6 6 a -\y + l-26) 



,-<f/2 



(87) 



quite similar to ( 1561) . While the relation a = b still holds true for the bcp and the bpcp 
with a < ac, that is no longer true for the bpcp at its multicritical point a = ac [5]. 



+ In the critical linear voter model [331 [92], the exact two-time autocorrelator [34] agrees in ID with 
the one of the ID Glauber-Ising model and for d > 2 with the one of the BCP. 
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field 


scaling dimension 


mass 






X 


M 







X 


-M 






x 2 


-2M 


T 


:=0 2 


it 


-M 


S 


:= 3 


Xt, 


-2M 


r 


= 3 2 


x r 


-M 



Table 3. Scaling dimensions and masses of some composite fields. 



We now review how the formalism of LSI can be adapted to these models, following 

®. 

4-5.1. Bosonic contact process From a field-theoretic point of view, the bcp can 
be described in terms of an order-parameter field and a response field 0, defined 
such that (0(t,r)) = (0(t, r)) = 0. The action [651 E3 [99] is again decomposed 
J\4>i 0] = i7o[0,0] + J'b[4 ) ,4 > } i^o a 'deterministic' part 



M,<f>] := dR du 



o 



(t>{2Md u - V 2 



and which is manifestly Galilei-invariant, whereas the 'noise' now contains two 
contributions and is described by 



^[0,0]:=-/// dR du 2 (0 + p o ) • (89) 

JR d JO 1 J 

where uncorrelated initial conditions are implied. 

In what follows, some composite fields will be needed. Together with their scaling 
dimensions and their masses, they are listed in table El For free fields one has x 2 = 2x, 
ix = 2x + x, = 3x + x and x^ = 3x + 2x, but not necessarily so for interacting fields. 
On the other hand, from the Bargman superselection rules we expect that the masses 
of the composite fields as given in table [3] should remain valid for interacting fields as 
well. 

As for magnets, the response function can reduced to a form which does not contain 
the noise explicitly [6]. The correlator becomes 

C(t,s;r -r') = (${t, r)0(s, r') exp ^-// J dR J du 2 (w, R)<f>{u, R) 



x exp I — fj,p I dR I du (J) 2 (R,u] 



o 



fipo JdRjdu (<f)(t, r)0(s, r')0 2 (w, R)^ (90) 
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+ y J dRdR' J dudu' (<j>(t, r)<p(s, r')T(u, R)T{u', R')) 

using the field T, see table [31 Hence the connected correlator is determined by three- 
and four-point functions of the noiseless theory. 

The noiseless three-point response can be found from its covariance under the ageing 
algebra aQt d [6] 

(<f>(t, r)<f){s, r')4> 2 {u, R)) = (t - sf'^ 2 (t - u)'^ 2 {s - u)'^ 2 

( Mir - R) 2 M(r'-R) 2 \ 
x exp M V_Jt) _ it) MuijVi)Q{t _ u)e(s _ u) (gi) 

with 

u [{s -u)(r-R)-(t- u)(r' - R)} 2 
Ul ~l (t-u)(s-u) 2 

u {(s-u){r-R)-{t-u){r'-R)] 2 

Vl ~S (t-uns-u) (92) 

and an undetermined scaling function Hence the critical bcp is described by a 
free field-theory, one can expect that x = x = d/2, hence X2 = d, xr = \d for 
the composite fields. Therefore the autocorrelator takes the general form C(t,s) = 
s 1_d/2 0i(t/s) + s 2 ~ d g2(t/s) and for d > 2 the second term, which comes from the four- 
point function, merely furnishes a finite-time correction to scaling. Now, if one chooses 
in eq. ([EE]) [87] 

tt 3 («i,Vi)=H (—--) (93) 

\Ui VxJ 

where H remains an arbitrary function, then 



d 



C(t,s)= -iip st +1 - x -^ 2 (t/s-l) 



jx 2 -x- 



,d 1 ; 



x / d6 [(t/s -6){l- ( t/S t J s 1 _ ] 26 ) (94) 



where the function \I/ is defined by 

Mw 



¥(w) = / dR exp — R 2 E(R 2 ) (95) 

and the LSI-prediction for the critical point (1551) is indeed recovered. The specific result 
for the bcp as listed in table [H follows for &(w) = $ w" 1 " n . 



4-5.2. Bosonic pair-contact process In this case, a new ingredient will be needed. The 
action now reads [65] J[4>, 4>] = Jo[<f>, 4>] + Jb[4>, 0] with the 'deterministic' part 



Joi 



dR 



dt 



<f>(2Md t - V 2 



(96) 
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Jb[<f>,4>] 



dR 



du 



2 j.3 



29 



(97) 



The representations of <XQt d or scf} d considered in sections 2 and 3 can no longer be used, 
since the equation of motion associated to Jq is non-linear, viz. 



2Md t (f>{t, r) = V 2 <P(t, r) - g4> 2 (t, r)0(t, r) 



(98) 



Rather, new representations of age d and sct) d must be constructed which take 
into account that g is a dimensionful quantity which transforms under local scale- 
transformations [96| [6] but the generators will not be written down explicitly here for 
the sake of brevity. It can then be shown that (i) eq. fUR) is Schro dinger- invariant for 
any valu^ of d and (ii) the Bargman supers election rules (3l\ ) still apply. 

Turning to the calculation of the autocorrelation function, one must consider five 
possible contributions. It turns out [6] that the results of tabled] for fc(y) in the bpcp 
can be reproduced from the single term 



C(t,s) 



ap 2 



dR / du (0(t,r)0(s,r)0 2 (w,i?) 



The required agt d - invariant three-point function now reads, where x = a + 1 
read off from the response function [6] 

'0(t,r)0(s, r')4> 2 (u,R) 



(99) 
d/2 is 



x exp 



(t- s) X S X2 (t-u) 



M(r- R? M (r 1 - R) 



5 X2 (s — u) z x ' 2 



u 



with 



m 
Pi 

«2 



2 t-u 2 

u l(s - u )(r-R)-(t- u)(r' - R)} 2 

1 (t - u)(s - u) 2 

u [(s — u ) (r - R) - (t- u)(r' - R)} 2 

s 



(t — u) 2 (s — u) 



a 



s 2 (t-u) 2 ' 

1 1 
t — u u 



02 



a 



i/.v 



s 2 (s — u) 2 '' 



03 



^ 3 (u 1 ,v 1 ,0 l ,0 2 ,0 3 )(lOO) 

(101) 
(102) 

(103) 
(104) 



a 1/y s 2 



where y is the scaling dimension of g. Consider the following choice for ^3 [6j 



^ 3 ( Ul , Vl , 0^02,03) = E 



(i 


-a 







03-V02~P~3 



(a-b) 



(105) 



where the scaling function 5 was already given in eq. (1931) . as for the bosonic contact 
process. We now have to distinguish the two different cases: first, if a < etc then a = b 

* If g is taken to be a dimensionless constant, it is a well-known mathematical fact that Schrodinger- 
invariance of (|98|) only holds true in 2D, see e.g. [40] . 
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[5] and we are back to the expressions found for the bcp. Second, if a = ac then 
a — b — d/2 — 1 for 2 < d < 4 and a — b = 1 for d > 4, respectively [5]. Identifying 
X2 = 2(6 — a) + <i, one finally obtains 



C(t,s) = s-\y- l)^)-"- 1 / d0[(t/s-0)(l 



la— 6 



X \I> 



t/s + 1- 29 
t/s- 1 



0(t/s-l) 



W/s-9){l- 



a—b 



(106) 



The specific results of table [T] for the bpcp are recovered if one uses the free-field form 



5. Growth models and an outlook towards z ^ 2 

A different class of non-equilibrium models considers the ballistic deposition of particles 
on a surface. The state of this surface may be described in terms of a height variable 
h(t,r). For irreversible deposition, the system clearly never arrives at an equilibrium 
state. Working in the frame co-moving with the mean surface height, the simplest 
kinetic equation one may write in the case without mass conservation is the well-known 
Edwards- Wilkinson (ew) model [35] 

d t h(t,r) = DV 2 h(t,r) +7](t,r) (107) 

However, if mass conservation must be taken into account, one might rather consider 
the Mullins- Herring (mh) model, see [801 H02j 

d t h(t, r) = -D(V 2 ) 2 h(t, r) + r](t, r) (108) 

Following [90], the following types of gaussian noise with vanishing first moment 
(rj(t, r)) = will be considered: 

(a) non-conserved, short-ranged (r](t, r)r](s, r')) = 2D5(r — r')5(t — s). 

(b) non-conserved, long-ranged (r](t,r)ri(s,r')) = 2D\r — r'\ 2p ~ d 5(t — s) and < p < 
d/2. 

(c) conserved, short-ranged (rj{t,r)rj(s,r')) = — 2DV 2 ,8{r — r')5(t — s). 
Then the following models were studied in [90] : 

(i) EWl: eq. (11071) with the non-conserved noise (a). 

(ii) ew2: eq. ( 11071) with the non-conserved, long-ranged noise (b). 

(iii) mhT: eq. f 1 1 8 p with the non-conserved noise (a). 

(iv) MH2: eq. fl 1 8 j) with the non-conserved, long-ranged noise (b). 

(v) MHc: eq. H108[) with the conserved noise (c) [9]. 
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To these, one may add the spherical model with a conserved order-parameter (model B 
dynamics) [TTJ [771 ESI E] which for d > 4 reduces to MHc. For 2 < d < 4 an explicit, 
if rather lengthy, expression for C(t,s) and its explanation in terms of LSI has been 
derived [9]. 

In the models EWl and MHc the noise is in agreement with detailed balance while for 
the other models it is not. Solving the linear equations H 1 71) and (11081) is straightforward 
and the two-time autocorrelations are listed in table [TJ The result quoted for the 
MHc model is only valid for d > 2 as stated; for d = 2 the scaling function becomes 
fc(y) — 2Z?ln[(y — l)/(y + 1)] [9j. Detailed simulations show that the correlation and 
response functions of the well-known Family model [37] and of a variant of it are perfectly 
described by the EWl model and hence should be in the same universality class [90] . 

5.1. Local scale-invariance for z = 4 

In line with the approach taken for z = 2, we shall consider the dynamical symmetries 
of the 'deterministic' part of these equations and then try to prove a generalization of 
the Bargman superselection rules such that the results can be extended to the full noisy 
system. Consider the 'Schrodinger equation' S&ip = where the 'Schrodinger operator' 
is defined as 

S, := -/j3t + ^d A r (109) 
lo 

The parameter /i is the analogue of the mass M. for the case z = 2 (and should not be 
confused with the exponent /i defined in section 3). Dynamical symmetries of this kind 
of simple linear equation have been constructed for any given z [50]. For the special 
case z = 4, some of the generators read in d = 1 dimensions as follows (see [90j [9] for 
the extension to d > 1) 

X_! = -d t 
Y-l/4 = - 9 r 

j £ 

X = - td t - -rd r - - 

X 1 = - t 2 d t - -trd r - |t - iir 2 d; 2 + 4 7 r<9- 3 - 6^d; 4 (110) 
F 3 /4 = - td r - 4/ird; 2 + 8 7 «9 r 7 3 

where 7 is a parameter and the 'derivative' 9" satisfies the formal properties <9" <9* = d^ +b 
and [d®,r] = ad^ 1 [50J. Then it is straightforward to check that 

[5 4 ,K.i/4] = [54, VI = [ s ^ Rii ' j) \ = > [Sa,X ] = -S 4 (111) 



and 

[S 4 ,X a ] = 



(112) 
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and the commutator of S4 with all other generators (IllOp vanishes. As before for 
Schrodinger and conformal invariance, this means the generators (11 101) act as dynamical 
symmetry operators provided that the scaling dimension of the solution if> of S^ip = 
is x = 3/2 — 27///. As before, a scaling operator is called quasiprimary [121 ES EH] 
if its infinitesimal transformation is given by ( 11 101) . Quasiprimary operators are now 
characterized by the triplett (x, 7///,//), rather than the pair (x,A4) as it is the case 
for Schrodinger-invariance. These 'quantum numbers' are connected to the critical 
exponents of the model at hand. It appears that x and 7//! are universal numbers, 
while /i (as is Ai) is dimensionful and must be non- universal. 

(2) 

The two-point function F built from two quasiprimary scaling operators and 
4> has the form [501 EH E] 

<V(t,r)0(s,r')) o = 5,,^7,7<W (t ~ s)~* /2 <p (|r - r'\{t - s^ 4 ) (113) 

where the scaling function is given by tp(u) = J2m=o c m^ m \u) where 

(-16/i) n r(l+m)rfn + f-^ 

y,M(u) - y^ 6 M M 4(n-l)+m+l ? fe M = V _ (n4) 

^ ^ r(4(n-l)+m + 2)r(f - j) 

Here, the coefficients c m are left arbitrary (although further constraints may be imposed 
on them by requiring that the scaling function <p(u) —>■ as u — > 00). 

In order to break time-translation invariance as required for applications to ageing, 
we introduce a time- dependent potential v(t) and consider the equation 



16// 



^-V* ( - V^> + v(t)<f^ + 77ft r) + j(t, r) (115) 



with a gaussian white noise i](t,r) and a perturbation j(t,r) for the calculation of 
responses. The potential can be eliminated through a gauge transformation [§] 



0(t,r) 



exp (-j^ J drv(r)Vl\ ^(t, r) (116) 



In contrast to the case z = 2, we now assume that for t — > 00 / *dru (r) ~ t r in order to 
ensure scaling behaviour for large times. This defines the new parameter F but we warn 
the reader that this parameter should not be confused with the one used in eq. 



Now, progress is possible if the deterministic part reduces to a linear equation. 
Then Wick's theorem holds which permits to derive analogous reduction formulae as in 
the case z — 2. For a fully disordered initial state with (4>(0, R)(f>(0, R)) = ao5(R — R) 
the response function is then explicitly given by [901 M 

R(t,s;r)=F^\t,s;r) 



3 



\( m ) 

,,k 

m=0 k,n=0 



(i-r /2 E^ m+1 EC 



Ko(t r — S 

16//(£ — s) 



h \ 



1/2 



k 



u 4(n-l)-2k (n7) 
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with the scaling variable u = \r\(t — s) -1 / 4 and the coefficients b ™^, which are given by 
h (m) _ 1 r(4(w-l)+m + 2) (m) 

° n ' k ~ k\T(2k-A(n-l)-m) n 1 J 

which is perfectly consistent with the available information from the several MH models 
introduced above (90] and also for the spherical model with a conserved order-parameter 

®- 

On the other hand, the correlation function becomes, with t > s 
C(t,s;r) = ^ [ di?F (3) (M,0;r,0,#) 

/ dr / dRF^ 3 \t,s,T;r,0,R) (119) 

./n JR d 



T 



16/i 

where the deterministic three-point function is 

F (£, s,r;r,r' , -R) := (0(t, r)0(s, r')0 2 (r, i?)) . If the deterministic part reduces to a 
linear equation, then Wick's theorem leads to [9J 

F (3) (t 1; t 2j t 3 ; n, r 2 , r 3 ) = 2F (2) t 3 ; n - r 3 )F (2) (t 2 , t 3 ; r 2 - r 3 ) (120) 

where F (2) can be taken from eq. (JHZJ. Carrying out the (rather long) calculation then 
allows to fully reproduce all results listed in tabled] for the models with z = 4. At least 
for the models studied so far with z = 4, the relation a = a' was empirically found to 
be satisfied. That is not surprising in view of the simple linear equation satisfied by the 
order-parameter field <fi. 

The MH models considered in this section are, together with the critical spherical 
model with a conserved order-parameter, the first analytically solved examples with 
z ^ 2 where local scale-invariance could be fully confirmed. These examples make it in 
particular clear that the height of the surface in growth processes is a natural candidate 
for being described by a quasiprimary scaling operator of local scale-invariance. 



5.2. 'Conformed invariance' for z ^ 1,2 ? 



We finish with a short speculation on how one might construct an analogy of the 
extension sch d — > conf d+2 discussed in section 2.3. By analogy, one might try to consider 
the 'mass' \x as a further variable and thus try to find analogues to the operators V- 
and N described in section 2.3. If we very naively generalise V- and iV from the case 
z = 2 and define 

V. = -id^ + 4ird t , N — —tdt — 1 — ^ (121) 

then it can be checked that both commutators [64, N] and [64, VI] vanish. This 
already implies that these two operators describe new, non-trivial symmetries and are 
appropriate candidates for a possible extension of the algebra freely generated from the 
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minimal set X , Xi, YL.1/4, ^3/4}- At this stage, the algebraic structure of such an 

extension is not clear at all; see [50] for the discussion of some other difficulties which 
may arise when one wishes to find the algebra satisfied by local scale-transformations 
with 2 ^ 1,2 and such that they can be symmetries of physically non-trivial equations 
- which in the context of S4 would mean that /1 7^ 0, 00. 

Of course, the central question whether an analogue of Bargman superselection 
rules exists is not even addressed. We hope to come back to this elsewhere, see 

6. Conclusions 

We have considered how local scale-invariance might be used in order to arrive at an 
explicit prediction for the two-time autocorrelation function of phase-ordering kinetics 
(as well as for non-equilibrium critical dynamics in those special cases where z = 2) 
capable to reproduce simulational data. In comparison with the two-time response 
function, whose functional form simply follows from the assumed covariance of the 
quasiprimary scaling operators of the ageing algebra age d , the calculation of C(t,s) 
does require conceptually important extensions. First, the explicit reduction from 
a stochastic Langevin equation to a deterministic equation is needed since a simple 
covariance argument for C(t,s) would have led to a vanishing autocorrelation because 
of the Bargman superselection rules, see [87]. Second, invariance under agt d alone does 
not seem to lead to any predictions for the autocorrelator, with the exception of the 
exponent relation Ac = \r, between the autocorrelation and the autoresponse exponents 
for the case of fully disordered initial conditions. Therefore, the third required extension 
concerns the dynamical symmetry algebra itself 

aQt d -> 5ci) d -> cotxf d+2 (122) 

which then allows to reduce the scaling function fc{y) to the linear combination of 
two known functions. We have also described some subleties in the precise definition 
of quasiprimary scaling operators and the several scaling dimensions (x,£, F) on which 
they depend. 

Besides several confirmations of the first two steps from analytically solved models 
with an underlying free-field theory, remarkably there is strong evidence from two- 
dimensional kinetic Ising and Potts models in favour of the last extension to a new type 
of conformal invariance, at least for quenches to T < T c , where z = 2. 

It might be useful here to list some of the important open problems. 

(i) A more systematic approach replacing the glueing procedure outlined in section 3.2 
(and the appendix) should be found, hopefully leading to full agreement between 
the theory and the data. 

(ii) The Galilei- and other local scale-invariances of non-linear partial differential 
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equations must be studied. Present results [951 E] require Galilei-invariance for 
all times which is too strong a requirement to be able to include the usual kind of 
Langevin equations [63] . A new kind of asymptotic Galilei-invariance, only valid in 
the scaling regime ([3]), would be enough and might be more flexible, but remains 
to be constructed. 

The importance of Galilei-invariance implies that models without a spatial structure 
such as the zero-range process [361 IE] cannot be expected to satisfy a local scale- 
invariance. On the other hand, it might be of interest to inquire whether driven 
diffusive systems [93] may possess some form of local scale-invariance. 

(iii) Analytical tests of local scale-invariance only exist for a rather simple kind of 
solvable systems, with an underlying free-field theory. The calculation of two-time 
observables in more rich integrable systems would certainly provide most useful 
information. 

(iv) Present tests of local scale-invariance have been limited mainly to systems with 
z = 2 or focused on the autoresponse whose form does not depend on directly on 
z. Therefore important aspects of the theory as outlined in [50] have not yet been 
tested. We hope that the consideration of phase-ordering in disordered Ising models 
[EU EE] and in certain models with long-range interactions j2U [TO] where z depends 
continuously on the control parameters of the model, will allow to do this. 

(v) A central point for the applicability to stochastic Langevin equations is the 
derivation of extensions of the Bargman superselection rules to z ^ 2. Work along 
these lines is in progress 

(vi) Is there a way to formulate an infinite-dimensional extension of the local 
scale-transformations considered here, in analogy with 2D conformal invariance 
at equilibrium [T2J [25]? At the time of writing, there merely exist some 
isolated mathematical results on equations with an infinite-dimensional dynamical 
symmetry, on the algebraic structure of the related Lie algebras and their vertex 
operator representations [291 Ell 1101] . What would be the physical signatures of 
such a symmetry which could be tested in specific models ? 

(vii) Dynamical scaling combined with reparametrization invariance has been considered 
for glasses where, because of logarithmic scaling, the effective dynamical exponent 
z — > oo, see [271 [28] and references therein. Will there be a way to relate these 
distinct approaches of local scaling ? 

All in all, local scale-invariance may offer an approach complementary to numerical 
simulation and the field-theoretical renormaliztion group (so far restricted to non- 
equilibrium critical dynamics [21]). The presently available evidence suggests the 
possibility of establishing a hidden dynamical symmetry which had not been suspected 
and, if proven correct, could simplify considerably the description of the non-equilibrium 
dynamics of many-body systems. 
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Appendix. 

We explain how to extend the description of the three-point function (<&(£)$(s)$ 2 («))o 
to the region where t ~ s and local scale-invariance can no longer be used. Throughout, 
we restrict to the case of phase-ordering (T < T c ) and we shall derive eq. f!76|) . Our 
physical criterion is the requirement that the two-time auto correlator be symmetric in 
t and s and furthermore be non-singular as t — s — > 0. From the from eq. (|49j) this leads 
to the following form of the scaling function \l/(p), since x — x 2 + 2£ — 2£ 2 = A^ — d/2 

oo 

y( P )=p X C- d / 2 J2*nP 2n (Al) 

71=0 

where the \l/ n are constants. Because of the Yeung-Rao-Desai inequality Ac > d/2 
[103] this expression will increase with p for sufficiently small values of p, whereas the 
LSI-prediction eq. ( 1721) gives a divergence at p = 0, unless a' < which in any case 
is incompatible with the explicit results in many models, see section 4 and [61, 62J. 
We expect (lAlj) to be valid for sufficiently small values of p and to change to the form 
predicted by LSI if p becomes larger than some minimal value e. For p sufficiently small, 
we may approximate (1721 as follows 

«-(,) * (i - 1 f - d d/ - 2 _ a > ,+ ■ ■ ■ 

+ ^M 2a ' +d / 2 - Ac + . . . (A2) 

which should be valid for p > e. Since in phase-ordering one usually treats systems of 
class S where a — a' = 1/z = 1/2 and A^ > d/2, all the terms retained here are needed 
to have a consistent expansion as p — > 0. 

If the cross-over point £ is sufficiently small, only the leading term in (1A1[) is needed. 
We construct a scaling function ^f(p) for all p > by requiring that the branches (]Alj) 
and flA2j) meet continuously at p = e. This condition fixes — ^o( £ ) an d we find, to 
this order 

T / \ , -in' i kj 2\q — d — 2a' i_2o' 
*o(e)=** 2a ~^ 1 + Ac _ rf/2 _ 2a / 

+ % p 1 M 2a ' +d l 2 - Xc e d/2 - Xc (A3) 
The required scaling function ^ is given by 



/■oo 

2 d ' 2 - 1 S d du«^- 1 e- jMwu *(w) 
2d/2_1 ^ (/°° dM Md/2_1 e_A1U,U ^ LSI ( M ) 

+ f dli M^ 2 " 1 (^o(^)« AC7 - d/2 - ^scal(w)) 
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=:2 d / 2 -% (J + h + I £ ) (A4) 

where Sd is the surface of the sphere in d dimensions, J ,i are proportional to ipo,i, ^lsi 
is the LSI-prediction given in (172j) and eqs. (1A2I) and (1A3I) were used. In the third line 
the range of integration was split into the intervals [0, e] and [e, oo] and the leading term 
from ( 1A1I) was used in the first subinterval. The terms 7o,i refer to the two independent 
solutions for in eq. ( 1721) . They can be calculated using the identity 

du u a ~ l e^ u 1 F 1 (A, B; -u) = 2 Fi \ A, a; B; ~j (A5) 

This already gives the first two lines in (1751) . The terms contained in I e are calculated 
from the identity 

i 

duu p - 1 e- uu = u- p 'y(P,u) (A6) 

where j(/3,u) is an incomplete gamma function p]. Combining the various terms and 
setting 

A := i> 2 d/2 - 1 S d M 2a '- Xc , B :=i; 1 2 d/2 - 1 S d M 2a '- Xc , E := Me (A7) 

we finally arrive at eq. (1761) in the text. 

Inserting this result into the two-time autocorrelation function eq. ( 15^1) . the 
amplitude = *(1) = lim^oo y x °/ 2 C(ys, s) is given by 

Coo = A {T(X C - 2a') 2 F 1 (2X C -d- 2a', X c - 2a'; 1 + X c - d/2 - 2a'; -1) 
+ E- 2a ' 1 (X c ,E)- 1 (X c -2a',E) 

+ E 1 - 2 *' - 2Xc ~ d ~„ 2a \ . [E^'-^Xc + 1 - 2a', E) - 7 (A C , E)] } 



1 + X c - d/2 - 2a' 
+ B {T{d/2) 2 F 1 (X C - d/2, d/2; 1 + 2a' + d/2 - X c ; -1) 

+E d/2-Xa l{Xcj E) _ 7(d/2j E) } (A8) 

which gives a useful technique to fix one of the free parameters A,B,E in terms of the 
other two. Indeed, in order to find empirically from numerical data, it is considerably 
more precise to plot C(ys, s)y~ x °^ 2 (y — l) Xc over against y, rather than the simplistic 
C(ys,s)y x o/ 2 . 

The opposite limit lim^! C(ys, s) = C(t, t) = C\ exists if 2a' < 1 and then becomes 
& = r(A c )2-^ Ue^' \l - ^ 1+ 2 ^~^~ 2 _ a/ 2Q/ (1 - cW)] + BE d / 2 - x o | (A9) 

which for finite E is indeed a finite constant. 

For a free field, one has Xc = d/2. If one takes A = and B > 0, one recovers 
#0) = BY{d/2)w' d / 2 , as it should be. 
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If one should find that the description of ^(p) for p < e merely by the lowest order- 
term is not sufficiently precise, further terms from flAlj) may be included (balanced by 
the higher-order terms of the expansion of \1/lsi(p)) and the coefficients \l/ n fixed by 
requiring also the continuity of the respective derivatives of ^(p) at p — e. 
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